Survival in equilibrium step fluctuations 
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We report the results of analytic and numerical investigations of the time scale of survival or 
non-zero-crossing probability S(t) in equilibrium step fluctuations described by Langevin equations 
appropriate for attachment/detachment and edge-diffusion limited kinetics. An exact relation be- 
tween long-time behaviors of the survival probability and the autocorrelation function is established 
and numerically verified. S(t) is shown to exhibit simple scaling behavior as a function of system size 
and sampling time. Our theoretical results are in agreement with those obtained from an analysis 
of experimental dynamical STM data on step fluctuations on Al/Si(lll) and Ag(lll) surfaces. 
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There has been much recent theoretical Q, and ex- 
perimental 0, interest in the persistence behavior of 
fluctuating steps on a vicinal surface. A persistence prob- 
ability P(t), defined as the probability that the position 
( "height" ) of the step edge at a point along a fluctuating 
step does not return to its initial value (at time t = 0) 
over time t is found 0, 0> m these studies to decay in 
time as a power-law, P(t) oc t~ e , for large t, where 9 is 
the so-called persistence exponent. Similar power-law be- 
havior of the persistence probability has also been found 
in experiments for other physical processes. It turns 
out that the precise definition of P(t) is absolutely crucial 
for the power-law behavior discussed in the recent sur- 
face fluctuations literature. If, instead of considering the 
probability of not returning to the initial position, one 
defines a survival probability, S(t), as the probability of 
the dynamical step height (at a fixed but arbitrary spatial 
location) not returning in time t to its average ( "equilib- 
rium") level, then, quite surprisingly, it was found in a 
recent experimental study [2| of thermal fluctuations of 
surface steps that S(t) actually manifests, in sharp con- 
trast to the power-law behavior of P(t), an exponential 
decay, S(t) oc exp(— t/r s ), at long times, where r s is the 
survival time scale. This exponential behavior of S(t) 
has remained theoretically unexplained. 

In this report we provide a definitive theoretical expla- 
nation for this exponential temporal behavior of the sur- 
face fluctuation survival probability using rigorous (an- 
alytical) arguments and direct (numerical) simulations. 
Survival and persistence turn out to be identical in prob- 
lems related to Ising spin dynamics, where one is inter- 
ested in the probability that a spin has not changed its 
sign (has not "flipped") up to time t This is due to 
the discrete nature of Ising spin dynamics, where a spin 
flip ensures a change of sign with respect to both the ini- 
tial and the average (or any other reference) value of the 
stochastic variable. In contrast, the continuous nature 



of surface fluctuation dynamics, where the step height 
is a continuous variable, leads to a fundamental qualita- 
tive difference between P(t) and S(t). This qualitative 
difference was noted as an experimental fact and an un- 
solved puzzle in Ref. Q for equilibrium step fluctuations 
on Al/Si(lll) @. We emphasize that the probabilities 
P(t) and S(t) provide completely different physical in- 
formation about surface step fluctuations: while persis- 
tence characterizes the universality class of the dynami- 
cal process through the persistence exponent 8, survival, 
as discussed in this paper, provides useful information 
about the physical mechanisms (and their characteristic 
time scales) underlying step fluctuations in the long-time 
limit. 

Another time scale that invariably enters experimental 
and numerical measurements of any statistical quantity 
is the sampling time 5t (the interval between successive 
measurements of the step position). An understanding 
of the effects of a finite St on the measured survival 
probabilities is necessary for comparing experimental and 
numerical results with theoretical predictions. 

In this paper, we present the results of a detailed 
study of the behavior of S(t) for two linear Langevin 
equations that describe 0, |3| step fluctuations under 
attachment / detachment ( "high-temperature" ) and edge- 
diffusion ("low temperature") limited kinetics. We first 
show analytically that if the equilibrium autocorrelation 
function C(t) of height fluctuations decays exponentially 
at long times, then S(t) must also decay exponentially 
with a time scale that is proportional to the correla- 
tion time (the time scale of the decay of C(t)). This 
prediction is verified from numerical simulations of the 
Langevin equations. The simulation results also provide 
information about the dependence of the measured sur- 
vival probability on the sampling time St. We show that 
the survival probability S(t, L, St) exhibits simple scaling 
behavior in both models. Finally, we use available exper- 
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physical systems that are believed to be described by 
these two Langevin equations and show that the experi- 
mental results are consistent with our predictions. 

High-temperature step fluctuations dominated by 
atomistic attachment and detachment at the step edge 
arc known 0, to be well described by the second-order 
non-conserved linear Langevin equation 



dh(x,t) _ T Q /3 d 2 h(x,t) 



Of 



k B T dx 2 



+ i](x,t). 



(1) 



Here, h(x,t) is the dynamical height fluctuation (po- 
sition of the step edge measured from its equilibrium 
value) at lateral point x along the step and time t, 
T a is the "step mobility", (3 is the step-edge stiffness, 
and n(x,t) is a nonconserved Gaussian noise satisfying 
{r](x, t)n(x' 7 1 1 )) = 2T a 5(x — x')5(t — t'). Low-temperature 
step fluctuations dominated by the step edge diffusion 
mechanism are, on the other hand, described by the 
fourth order conserved Langevin equation 



dh(x,t) T h f3d i h{x,t) 



dt 



k B T dx 4 



r) c (x,t), 



(2) 



with ( Vc {x, t)r) c (x', t')) = -2T h V 2 J(x - x')S(t - t'). 

Space- and timc-dcpcndcnt correlation functions of 
height fluctuations in these two linear equations may be 
calculated 0,0 easily by Fourier transforms. We assume 
that height fluctuations are measured from the spatial 
average of h(x,t), so that the k = Fourier component 
of h{x,t) is zero at all times. The autocorrelation func- 
tion of h(k,t), the Fourier transform of h(x,t), has the 
following form in the long-time equilibrium state: 



(h(k,h)h{-k,t 2 )) 



k B T 

f3k 2 



exp(-r^ z |i 1 -t 2 \/k B T), 



(3) 

where z — 2, T = T a for Eq.JTJ, and z = 4, T = T h for 
Eq.lJSJ. The autocorrelation function of height fluctua- 
tions at equilibrium is then given by 



C{t) = (h(x,h)h{x,t 2 )) 

2k B T f 00 dkexp(-T/3k z \t\/k B T) 



(4) 



where t = t\ — t 2 and k m \ n = 2tt/L for a finite system 
of linear dimension L. This implies that C(t) exhibits 
an exponential decay at long times, C(t) oc exp(— t/r c ), 
where the correlation time r c is equal to k B TL z / (2tt) z TP. 

There are other physical mechanisms that may lead to 
an exponential decay of C(t) at long times. The fluctu- 
ations of a particular step arc affected by its interaction 
with neighboring steps. These interaction effects, which 
are negligible at relatively short time scales if the spacing 
between neighboring steps is large, may become impor- 
tant at long times. If one assumes that the step fluctuates 



in a harmonic confining potential yj, (A/2) J h 2 {x)dx, 
then one obtains an additional term, —XT a h(x,t)/k B T, 
in the right-hand of Eq.JQ, and \T h V 2 x h{x, t)/k B T in 
Eq.(j2J. The function C{t) is then given by 



C(t) = 2k B T 



dk exp[-r(/3fc 2 + \k z - 2 )\t\/k B T)} 



2tt 



(ik 2 + A 



(5) 

This, again, leads to an exponential decay of C(t) at long 
times, with the correlation time r c a function of A and L. 

There exists a rigorous theorem 0] that states that if 
the autocorrelation function C(t) of a stationary Gaus- 
sian process decays exponentially in time, then its sur- 
vival probability S(t) must also decay exponentially for 
large t, S(t) oc exp(— t/r s ), with the survival time scale t s 
proportional to the correlation time r c . The constant of 
proportionality, c = t s /t c , which must be less than unity 
and independent of the system size L, is usually nontriv- 
ial, being determined by the full functional form of C(i). 
Since the height fluctuations in our models represent a 
stationary Gaussian process at equilibrium, this rigorous 
result applies for the survival probability of these fluctu- 
ations. Thus, we arrive at a very general, exact result 
that the survival probability of equilibrium step fluctu- 
ations should decay exponentially at long times if the 
autocorrelation function does so. This is the first impor- 
tant result of this article. Further, measurements of the 
ratio c of the two time scales may provide valuable infor- 
mation about the nature of the processes involved in the 
step fluctuations. 

We have investigated these aspects in a detailed nu- 
merical study in which a simple Euler scheme is used 
to numerically integrate spatially discretized versions of 
Eqs.QI and (0. All the results reported here were ob- 
tained in the equilibrium regime. S(t) is measured as the 
probability that the height fluctuation hi at a particular 
site i does not cross the average step height (which is 
conveniently chosen as the "zero" of the height stochas- 
tic variable) over time t, averaged over all sites and many 
(~ 10 4 — 10 6 ) independent runs. C(t) is calculated ex- 
actly using discretized versions of Eqs.Q and 

Typical results for C(t) and S(t) are shown in FigQ] 
(a-c) for Eq.JTJ and in Fig0 (a-c) for Eq.© [A = in 
both cases]. As indicated in the figures, we used several 
different values of the sampling time St in the measure- 
ment of S(t) [C(t) is, of course, independent of St]. It 
is clear from the plots that both C(t) and S(t) decay 
exponentially at long times. The time scales r c and r s 
are extracted from exponential fits shown in the semi-log 
plots as dashed straight lines. The dependence of t c on 
L is given exactly by t c (L) = (L/2tt) z . The results for 
different values of L shown in Figs. and indicate that 
t s also increases rapidly as L is increased. However, we 
find that the calculated values of r s extracted from the 
S(t) data obtained for different L using the same sam- 
pling time St exhibit small but clear deviations from the 
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FIG. 1: S(t) and C(t) for the Langevin equation of Eq.(l). 
The dashed lines are fits of the long-time data to an expo- 
nential form. In panels (a-c), the uppermost plots show the 
data for C(t). Panel (a): L = 100, St = 0.625. Panel (b): 
L = 200, St = 2.5. Panel (c): L = 400, St = 10.0 (upper plot) 
and St — 1.0 (lower plot). Panel (d): Finite-size scaling of 
S(t, L, St). Results for S for 3 different sample sizes with the 
same value of 8t/L z (z = 2) are plotted versus t/L z . 



FIG. 2: S(t) and C(t) for the Langevin equation of Eq.(2). 
In panels (a-c), the uppermost plots show the results for C(t). 
The dashed lines are fits of the long-time data to an exponen- 
tial form. Panel (a): L = 20, St = 1. Panel (b): L = 40, 
St = 16 (upper plot) and L — 40, St — 1 (lower plot). Panel 
(c): L — 80, St — 256. Panel (d): Finite-size scaling of 
S(t, L, St). Results for S for 3 different sample sizes with the 
same value of St/L z (z — 4) are plotted versus t/L z . 



expected proportionality to L z . These small deviations 
result from a weak dependence of S(t) on the sampling 
time St. As shown in FigQJ; and Fig|2), the rate of the 
exponential decay of S(t) at large t depends weakly on 
the value of St used in the measurement of S. This is in 
accordance with the analytic predictions of Ref . pj . Since 
the only time scale in the problem is r c (as mentioned 
above, r s should be proportional to r c ), the dependence 
of S(t) on the sampling time St should involve the scal- 
ing combination St/r c . Since t c (L) oc L z in our models, 
this argument suggests that the sampling time should be 
chosen to be proportional to L z if the survival probabil- 
ities for different values of L are to be tested for scaling. 
Indeed, as shown in Fig|TJi, the values of S(t) obtained 
for L = 100, 200 and 400 using St = 0.625, 2.5 and 10.0, 
respectively, (so that St oc L z with 2 = 2) all fall on the 
same scaling curve when plotted as functions oit/L z with 
z = 2. As shown in Fig|2i, a similar scaling collapse is 
obtained for Eq.Q. Here, the sampling times for differ- 
ent L are chosen to be proportional to L 4 , and the best 
scaling collapse is obtained when the data for different 
L are plotted against t/L z with z ~ 3.95. These results 
establish that the full function Sit, L, St) (not just the 
asymptotic long-time part) has the scaling form 

S(t,L,St) = f(t/L z ,St/L z ), (6) 

where the function f(x,y) decays exponentially for large 
values of x and the rate of this decay increases slowly as y 
is decreased. This finite-size scaling behavior of S, which 



represents the second important result of our study, is 
similar to that found 01 f° r the persistence probability 
in a coarsening system. However, the dependence on the 
sampling time, essential in our scaling considerations, 
was not analyzed in Ref.|l2|. 

We have also studied the behavior of C(t) and S(t) 
for Eq.Q when the value of r c is primarily determined 
by the presence of a nonzero A associated with step-step 
interaction (cf. Eq. ||SJ above). By varying A and St 
for a system with L — 400, we find that S(t,X,8t) ex- 
hibits excellent scaling behavior as a function of t/r c if 
the quantity St/r c is held constant. Therefore, we con- 
clude that S is a function of the scaling variables t/r c 
and St/r c , irrespective of the origin of the finite value of 
the correlation time r c . 

For A = 0, the ratio c = r s /r c for Eq.Q decreases from 
about 0.57 to about 0.41 as the ratio St/r c is decreased 
from 0.025 to 2.5 x 10~ 4 , indicating that c ~ 0.4 in the 
St —> limit. For relatively large L and nonzero A where 
t c is determined primarily by the value of A, we find that 
c ~ 0.47 for 8t/r c — 0.025. The difference between the 
values of c for the same value of St/r c in the two cases 
reflects the expected dependence of c on the details of 
C(t). For Eq.© with A = 0, the value of c decreases from 
about 0.44 to about 0.30 as St/r c is decreased from 0.01 
to 6 x 10~ 4 . The qualitative behavior of c as a function 
of St/r c is similar in all the cases we have considered, and 
is consistent with the general predictions of Ref. • 

We have also used dynamical STM data to calculate 
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FIG. 3: S(t) and C{i) for two experimental systems. The 
dashed lines are fits of the long-time data to an exponential 
form. Panel (a): Al/Si(lll) at T = 970K. Panel (b): Ag(lll) 
at T = 450K. 

C(t) and Sit) for two experimental systems: Al/Si(lll) 
at relatively high temperatures, which is believed |2j, 
Eii Q| to provide a physical realization of Eq.JU, and 
Ag(lll) at relatively low temperatures where the step 
fluctuations are expected 0, [ljj to be governed by the 
conserved Eq.JSJ). Some of the results of this analysis are 
shown in FigO For Al/Si(lll) at 970K we find exponen- 
tial decay of both C(t) and S(t). The value of the ratio 
c obtained from the estimates of t c and t s is close to 0.5. 
This value is right in the middle of the range of values 
of c obtained from our numerical study of Eq.Q. The 
Ag(lll) at 450K data is characterized by c ~ 0.34, which 
is again in the range of values obtained in the numerical 
study of Eq.0. We, therefore, conclude that the avail- 
able experimental data on S(t) and C(t) are consistent 
with our theoretical results. 

Experimental data on the Al/Si(lll) system are avail- 
able at several temperatures between 770K and 1020K. 
As reported in Ref.|2|, Sit) decays exponentially at long 
times at all these temperatures, with t s decreasing from 
3.6s to 0.9s as T is increased from 770K to 1020K. Using 
these values of t s (actually, the corresponding values of 
t c obtained from the relation r s /r c ~ 0.5) together with 
the values of the parameters V and [3 obtained from other 
measurements Qj3, , we have calculated an "effective 
length" L e f / that would lead to the observed finite value 
of r c if it resulted from a finite length of the sample. The 
value of L e ff is found to decrease from 4020A to 389A 
as T is decreased from 1020K to 770K. These values are 
much smaller than the nominal step lengths in the ex- 
perimental sample. The observed T-dependence of L e f f 
is inconsistent with the possibility that the finite values 
of t c are due to a nonzero value of the parameter A: the 
length scale associated with A should increase || as T 
is decreased. It is possible that L e // is a measure of 
the typical length of a step edge between adjacent points 
that are held fixed by some kind of pinning centers. Since 
pinning becomes more effective at low T, this mechanism 
would provide a qualitative explanation of why L e f f de- 
creases as T is reduced. Yet another possibility is that 



L e f f is a measure of the length scale over which step edge 
fluctuations are effectively equilibrated. 

To conclude, we have shown analytically and numer- 
ically that the survival probability of equilibrium step 
fluctuations on vicinal surfaces decays exponentially at 
long times, and have established a relation between the 
time scales characterizing the exponential decay of the 
survival probability and the autocorrelation function. 
Our theory explains the puzzling experimental finding 
of an exponential decay of S(t) reported in Ref. 0. We 
have also shown that the survival probability exhibits 
simple scaling as a function of the system size and the 
sampling time, which plays a very important role in the 
measurement of Sit). 
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